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Abstract. Observational appearance of a neutron star in the subsonic propeller state which is a companion of 
a wind-fed mass-exchange close binary system is discussed. During the subsonic propeller state (which was first 
introduced by Davies et al. 1979) the neutron star magnetosphere is surrounded by a spherical quasi-static plasma 



envelope, which is extended from the magnetospheric boundary up to the star accretion radius. The energy input 
to the envelope due to the propeller action by the neutron star exceeds the radiative losses and the plasma 
temperature in the envelope is of the order of the free-fall temperature. Under this condition the magnetospheric 
boundary is interchange stable. Nevertheless, I find that the rate of plasma penetration from the envelope into the 
magnetic field of the neutron star due to diffusion and magnetic field line reconnection processes is large enough 
for the accretion power to dominate the spindown power. I show that the accretion luminosity of the neutron 
star in the subsonic propeller state is L a ~ 5 x 10 30 -j- 10 33 Mis ergs -1 , where M15 is the strength of the normal 
companion stellar wind which is parametrized in terms of the maximum possible mass accretion rate onto the 
neutron star magnetosphere. On this basis I suggest that neutron stars in the subsonic propeller state are expected 
to be observed as low luminous accretion-powered pulsars. The magnetospheric radius of the neutron star in this 
state is determined by the strength of the stellar wind, M c , while the accretion luminosity is determined by the 
rate of plasma penetration into the star magnetosphere, M a , which is M a <C M c . That is why the classification 
of the neutron star state in these objects using the steady accretion model (i.e. setting M a = M c ) can lead to a 
mistaken conclusion. 
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1. Introduction 



As was first recognized by Shvartsman (1970), three states 
of a rotating, magnetized neutron star in a close binary 
system can be distinguished: ejector, propeller and accre- 
tor. This classification reflects three different evolutionary 
stages and three different mechanisms of energy release 
responsible for the neutron star emission. The sequence 
of states which a magnetized neutron star in a wind-fed 
mass-exchange binary system follows as it spins down from 
the initially very short periods can be expressed in the 
form of the following chain: ejector — > propeller — > accre- 
tor. 

Neutron stars in the state of ejector are known as spin- 
powered pulsars, in which the spindown power, 
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dominates the star energy budget and is spent predom- 
inantly on the generation of magneto-dipole waves and 
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particle acceleration. Here fi is the neutron star magnetic 
moment, (3 is the angle between the rotational and mag- 
netic axes, lu — 2tt/P s is the star angular velocity, P s is 
the neutron star spin period and c is the speed of light. 

The state of a neutron star is classified as accretor if 
the star luminosity is dominated by the accretion power, 

L a = M a ^. (2) 

Here M ns and R ns are the mass and the radius of the 
neutron star, respectively, and M a is the mass accretion 
rate onto the neutron star surface. 

The intermediate state is called propeller. The neces- 
sity to introduce this state can be justified as follows. The 
pulsar-like spindown ceases when the pressure of the ma- 
terial being ejected by the neutron star, 
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can no longer balance the pressure of the surrounding gas, 
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Here poa is the density of the surrounding gas, V Te \ is the 
relative velocity between the neutron star and the sur- 
rounding gas, and R a is the accretion radius of the neu- 
tron star, 



R a — 



2GM n 
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Combining Eqs. (|) and (@) one finds that this occurs 
when the neutron star spin period P s — P mc \, where 



P md =0.3 ^V" 1/4 (sin/3) 1 / 2 
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Here /x 30 = ^/10 30 Gcm 3 , V 8 = V rcl /10 8 cms" 1 and M c is 
the strength of the normal companion stellar wind which is 
parametrized in terms of the maximum possible accretion 
rate following^] Davies et al. ( L979 ) 



rel • 



(7) 



On the other hand, for the steady accretion onto the 
neutron star surface to realize the canonical magneto- 
spheric radius, which is defined by equating the ram pres- 
sure of the inflowing material with the magnetic pressure 
due to the dipole field of the neutron star, 



R m = 



2/7 



should be smaller than the corotational radiusQ, 
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Otherwise the neutron star proves to be in the centrifugal 
inhibition regime (i.e. the centrifugal acceleration at the 
magnetospheric boundary, u> 2 R m , dominates the gravita- 
tional acceleration, GM ns /Rf n ) and hence, the plasma ac- 
cretion onto the star surface is impossible. Putting Eq. (|J) 
to ([)[) one finds that the condition r m ^, R COI is satisfied 
only if M c M ca , where 
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Here m is the mass of the neutron star expressed in units 
of Mq. However so huge a mass transfer rate can never 
be realized in a wind-fed mass-exchange close binary sys- 
tem (Mea exceeds the maximum possible rate of mass loss 



1 There should be no confusion between M c , which is the 
mass of the surrounding material interacting with the neutron 
star moving through the wind in a time unit and M a , which 
is the mass accretion rate onto its surface. In the general case 
M c ^ Ma,. The case of M c = Af a corresponds to the steady 
accretion process, when all material captured by the neutron 
star reaches the star surface. In this particular case M c is called 
the r ate of mass capture by the neutron star (Bondi & Hoyle 
1944). However, as it will be shown below, it has a different 



physical meaning if the neutron star is in the state of propeller 
2 This is a necessary but not a sufficient condition 



by O-type supergiants by more than five orders of mag- 
nitude !). That is why the direct ejector — ► accretor state 
transition under the conditions of interest cannot be real- 
ized and the evolutionary track of the initially fast rotat- 
ing magnetized neutron star must contain an additional 
spindown state which is traditionally called propeller. 

A neutron star in the state of propeller is spinning 
down due to the interaction between its fast rotating mag- 
netosphere and the surrounding material. The mechanism 
of this interaction has been investigated in detail by Davies 
et al. ( 1979] ) and Davies & Pringle (1981) who have shown 
that during the propeller state the star magnetosphere is 
surrounded by a spherical quasi-static envelope in which 
the plasma temperature is of the order of the free-fall tem- 
perature, 



T{R)~T S (R) 
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Here m p is the proton mass and k is the Boltzmann con- 
stant. The rotational energy loss by the neutron star is 
convected up through the envelope by the turbulent mo- 
tions and lost through its outer boundary. The structure 
of the envelope and the spindown rate of the star depend 
on the value of the ratio: 

uR m 



V s (R m )' 



where i? m is the magnetospheric radius of the neutron star 
and V s (R m ) is the sound speed in the envelope plasma, 
which according to Eq. ( |i"l| ) is of the order of the free-fall 
velocity, Vg: 



V s {R m ) ~ V B (R„ 
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On this basis, Davies et al. ( 1979| ) distinguished three 
sub-states of the propeller statey 1 : (b) very rapid rotator 
(k^> 1); (c) supersonic propeller (k 1) and (d) subsonic 
propeller (n < 1). 

In cases 'b' and 'c' the magnetospheric radius of the 
neutron star exceeds its corotational radius. The neutron 
star is in the centrifugal inhibition regime and no accre- 
tion onto its surface occurs. The spindown power domi- 
nates the energy budget and is comparable to the radiative 
losses of the envelope. 

When the neutron star spin period decreases below the 
critical period, 
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the corotational radius reaches the magnetospheric radius. 
Under the condition P s > P cc \ the wind plasma, penetrat- 
ing from the base of the envelope into the star magnetic 
field, is able to flow along the magnetic field lines and to 
accrete onto the star surface. However the rate of plasma 



According to their classification, case 'a' corresponds to the 
pulsar-like spindown, i.e. the ejector state. 
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penetration into the star magnetic field remains essentially 
smaller than M c if the magnetospheric boundary is sta- 
ble with respect to interchange instabilities (e.g. Rayleigh- 
Taylor instability). According to Arons & Lea (1976a) and 
Eisner & Lamb (1976), the onset condition for the inter- 
change instability of the magnetospheric boundary reads 



T(P m ) < T cr ~ 0.3 T ff . 



(14) 



This means that the steady accretion process under the 
condition P s > P c d can be realized only if the cooling of 
the envelope plasma (due to the radiation and convective 
motions) dominates the energy input to the envelope due 
to the propeller action by the neutron star. 

Inv estigating this particular state, Davies & Pringle 
( |l98l| , hereafter DP81) have shown that the cooling pro- 
cesses in the envelope are effective if the spin period of the 
star exceeds a so-ca lled break period Pb r , which according 
to Ikhsanov (2001a) is 



Pbr ^ 450 fil 6 /21 M 



r -5/7 
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DP81 have assumed that during the epoch when P c( j < 
P s < Pbr, no accretion onto the star surface is realized 
(i.e. M a = 0) and the star energy budget is dominated by 
the spindown power, which in the case of this - subsonic 
propeller - state is 



L ssp = 6.4 x lO^egrs- 1 /xf m^Pr 3 



25 ' 



(16) 



where P25 = P s /25s. 

At the same time, the assumption M a = under the 
condition P m < P cor is not obvious since the interchange 
instability is not the only mode by which the accreting 
plasma can enter the neutron star magnetic field. In this 
paper I show that the accretion luminosity of the neutron 
star in the subsonic propeller state can be comparable or 
even exceeds the spindown power if the plasma penetra- 
tion into the magnetosphere is governed by the magnetic 
field lines reconnection process (section 2). In this case the 
star can be observed as a low luminous accretion-powered 
X-ray pulsar rather than the spin-powered propeller. The 
results are discussed and summarized in section 3. 



Eq. |l6|) dominates the radiative losses from the envelope. 
Under this condition the envelope can be considered as a 
static adiabatic atmosphere in which the plasma pressure 
p oc P~ 5 / 2 , the sound speed c s cx P -1 / 2 and the num- 
ber density n cx P~ 3 / 2 . The plasma pressure at the outer 
edge of the envelope is equal to the ram pressure of the 
surrounding gas. The formation of the envelope in the first 
approximation prevents the stellar wind from penetrating 
into the accretion lobe of the neutron star (i.e. the sphere 
with the radius R a ). In this case stellar wind overflows 
the outer edge with the rate M c compressing and heating 
the envelope plasma (for discussion see DP81). 

The equilibrium shape of a neutron star magneto- 
sphere surrounded by a hot hydrogen envelope has been 
reconstructed by Arons & Lea (1976a), Eisner &Lamb 



(1977) and Michel (1977). All these authors have argued 
that the magnetosphere in the considered case tends to 
be closed and its boundary is convex towards the accret- 
ing plasma, with two cusp points situated on the mag- 
netic axis of the neutron star. It is rotationally symmetric 
about the magnetic axis and reflection symmetric about 
the equatorial plane. 



2.1. The rate of plasma entry into the star 
magnetosphere 

Formation of the magnetosphere, to a first approximation, 
prevents the envelope plasma from reaching the neutron 
star itself. That is why the rate of mass accretion onto 
the neutron star surface strongly depends on the rate of 
plasma penetration into the magnetosphere at its bound- 
ary. As was been mentioned in the introduction, under the 
condition P s < Pb r , the cooling of the envelope plasma 
is not effective and hence the magnetospheric boundary 
is interchange stable. Under this condition the modes by 
which the envelope plasma can enter the star magneto- 
sphere are (i) turbulent diffusion and (ii) reconnection of 
the magnetic field lines. 



2. Accretion luminosity of a neutron star in 
'subsonic propeller' state 

I consider a close binary system in which a rotating, mag- 
netized neutron star is orbiting around a high-mass main 
sequence companion which underfills its Roche lobe and 
loses mass in the form of stellar wind. The neutron star, 
which moves through the wind of the normal companion, 
is assumed to be in the state of subsonic propeller. This 
means that the neutron star magnetosphere is surrounded 
by a hot plasma envelope (T « Tg) which is extended from 
the magnetospheric boundary, P m ^ P C or, up to the ac- 
cretion radius of the neutron star, P Q . The spin period of 
the neutron star lies within the interval P c d < P s < Pb r) 
and thus, the rate of energy input to the envelope due 
to the propeller action by the neutron star (expressed by 



2.1.1. Diffusion 

The conductivity of the envelope plasma can be approxi- 
mated by that of a fully ionized gas: a ~ 10 7 T e 3/2 , where 
T e is the electron temperature. In our case, the value of 
conductivity is very high but, nevertheless, finite (a < 00). 
That is why a diffusion skin layer (magnetopause) at the 
magnetospheric boundary is formed. The thickness of the 
magnetopause in the general case depends on the velocity 
and the time of plasma diffusion across the magnetic field 
lines, but in any case it is larger than the Larmor radius of 
ions. Under the condition P cc j J> P s , the gravitation force 
applied to the plasma in the diffusion layer dominates the 
centrifugal force. Therefore the plasma penetrating the 
layer flows along the field lines toward the magnetic poles 
of the neutron star and accretes onto the star surface. 
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The rate of plasma flow into the magnetopause from 
the envelope due to diffusion is 



Ma 



47Ti4 p(R m ) Vdiff- 



(17) 



p(Rm) is the plasma density at the base of the envelope 
just over the magnetosphere boundary. It can be evaluated 
equating the plasma ram pressure and the magnetic field 
pressure at the boundary: 



p(Rm) 



AnGM DS R^ 
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Vdiff is the diffusion velocity which, according to Ikhsanov 
& Pustil'nik (1996), can be expressed in the following form 



Vdiff — v Asff/iff- 



(19) 



Here D c h is the effective diffusion coefficient and tff — 
Rm 2 /VGM ns is the free-fall time, which is the character- 
istic time of the plasma flow from the layer to the star 
surface. 

In the general case the value of D e g is obviously limited 
by the following relation: 

D CI <>D cS <,D B , 

where D cr = c 2 /47rcr and Z?b = C^T(i) r h(i) are the coef- 
ficients of the classical and Bohm diffusion, respectively, 
Vr(i) and rvn are the thermal velocity and Larmor radius 
of ions, and ( is the diffusion efficiency. 

The minimum rate of plasma entry into the magneto- 
sphere can be obtained putting D c g ~ D c \. In this case 
one finds 



M cl ~ 10 6 p% 7 m- 15 / 14 
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The maximum inflow rate is realized if the diffusion 
process is governed by drift-dissipative instabilities (i.e. 
Bohm diffusion) 



(ckTj(R n 
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where Tj is the ion plasma temperature at the base of the 
envelope. Combining Eqs. (|l7]- 21) yields 




n V omlO -1 A/ 2 -1/14 1/7 

M B - 310 1 gs C0.1 M30 m 



where Co.i — C/0-1- This normalization of the diffusion 
efficiency is adopted following the results of experiments 
on the nuclear fusion (e.g. Hamasaki et al. 1974) and the 
measurements of the solar wind penetrating the magneto- 



sphere of the Earth (Gosling et al. 1991). 

Hence, the total rate of mass flow into the magneto- 
sphere due to the diffusion process is Md = M c \ + Mb ■ 



2.1.2. Reconnections 



According to Eisner & Lamb (1984) the rate of plasma 
entry into the magnetosphere due to reconnection of the 
magnetic field lines can be expressed as 



GM ns i? m AttRI 



V a , 



(23) 



where a r is the efficiency of the reconnection process, A r = 
47ri2 m A m is the effective area of the reconnection region 
and V a = B/^/Aup is the Alfven velocity, which in the 
boundary layer at the base of the hot envelope is of the 
order of Vs(R m ). 

The values of a r and A r depend on the current sheet 
parameters and the magnetic field configuration in the 
accretion flow. Investigations of reconnection processes in 
solar flares and in the Earth's magnetopause su ggest the 
average value of a T ~ 0.1 (see Priest & Forbes [2000 and 
references therein). 

The effective area of the reconnection region depends 
of the scale of the magnetic field in the accreting plasma 
over the magnetospheric boundary. It is limited by the 
initial inhomogeneity of the accretion flow and/or by the 
fragmentation of plasma over the rotating magnetospheric 
boundary due to the Kelvin- Hchnholtz instability. In the 
latter case the average scale of plasma vortic es of t he em- 
bedded field is A m ~ 0.1i? m (see Arons & Lea 1976b; Wang 
& Robertson 1985). On the other hand, the scale of the 



plasma inhomogeneity at the base of the envelope can be 
evaluated by the size of the convective cell (i.e. by the 
height of the homogeneous atmosphere over the magneto- 
spheric boundary, h g ). In our case (T <; 0.3Tff) one finds 
h g ^ 0.3i? m and, correspondingly, the ratio of the area of 
a convective cell to the total area of the magnetospheric 
boundary is (h g /R m ) 2 ~ 0.1. In this context the estimate 
A m ~ 0.1i? m seems to be appropriate. 

Under these conditions the rate of plasma penetration 
into the neutron star magnetic field due to the magnetic 
lines reconnection reads 



M n 
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Thus, the rate of plasma penetration into the magneto- 
sphere of the neutron star in the state of subsonic propeller 
can be expressed as 



M a = M d + M T , 



2.2. Accretion luminosity 



(25) 



Under the condition R m < R COI plasma penetrating the 
magnetic field of the neutron star is accreted onto its sur- 
face. This indicates that the observed total luminosity of 
the neutron star in the state of subsonic propeller is 
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bol 
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The lower limit to the accretion luminosity, L a , can be 
estimated assuming the plasma entry in the stellar mag- 
netic field from the base of the envelope due to diffusion 
(see Eq. ||): 



I 30 ^V. 2 ,," 1 / 14 ™ 8 / 7 Rr 1 M}}/ u erg s" 1 , (27) 

where M15 is the strength of the stellar wind expressed in 
units of 10 15 gs _1 . 

Combining Eqs. © and © one finds that the ac- 
cretion power L a ,diff dominates the star energy budget if 
is ^ -Pdda, where 



P 



dda 



63 #« ^/ 42 m- 5 /^V 3 Mr 5 11/42 



(28) 



If the plasma penetration the magnetosphere of a neu- 
tron star is governed by the reconnection of the magnetic 
field lines, the accretion luminosity of the star is 



10 



:i:S 



out m Rq 1 A .iMi 5 ergs l , 



(29) 



where ao.i = Qi r /0.1 and A0.1 = A m /0.1i?. m . In this 
case the accretion power already dominates the spindown 
power at P c d if M c > M r d a , where 

M rda ~ 6 10 13 /4 ™~ 2 i?6 aai A o.i ^ gs -1 - (30) 

In both cases the neutron star will be recognized as a 
low luminous accretion-powered pulsar in which the radius 
and the temperature of polar caps areP| 

1/7 
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amount of material coming into the envelope through its 
outer boundary. This means that the radial drift of plasma 
through the envelope with the rate M a is expected. Hence 
the question about the structure of this modified envelope 
arises. 

The characteristic time of the accretion process in the 
envelope can be estimated as 



td T — R/Vdr, 



(31) 



where Vd r is the average velocity of plasma radial 
drift through the envelope, which under the condition 
M a = const, can be evaluated as 



Ma 

Mr 



■Vs. 



(32) 



On the other hand, the time scale on which the incoming 
energy is redistributed within the envelope (i.e. the time of 
envelope relaxation to the equilibrium state) is determined 
by the characteristic turbulence time 



t t = R/V t . 



(33) 



Hence the quasi-static approximation can be applied to 
the modified model of the envelope if tt <C tdr- Combining 
Eqs. ( |3l| ) and ( |33|) I find this condition to be satisfied if 
the spin period of the star is P s -C P qs , where 



Pa, 
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This indicates that under the conditions of interest the 
structure of the modified envelope (with the radial plasma 
drift) can be considered within the quasi-static approxi- 
mation and thus, the model presented by DP81 can be 
applied. 



2.3. Quasi-static envelope approximation 

The model of the envelope which is surrounding the mag- 
netosphere of a neutron star in the subsonic propeller state 
has been constructed by DP81 under the assumption that 
no accretion onto the star surface occurs. In this case the 
envelope can be considered as quasi-static, i.e. the average 
rate of mass transfer through the envelope in the radial 
direction is assumed to be zero. Under this condition the 
plasma of the stellar wind does not penetrate into the ac- 
cretion lobe of the neutron star but overflows the outer 
edge of the envelope with the rate M c , compressing the 
envelope plasma. 

One faces however a slightly different situation taking 
into account that the "magnetic gates" during the sub- 
sonic propeller state of the neutron star are not closed 
completely As has been shown above, the rate of plasma 
inflow into the interchange stable magnetosphere due to 
diffusion and the magnetic field line reconnection is M a ^ 
0. In this case the envelope remains in an equilibrium state 
if the amount of material inflowing from the base of the en- 
velope to the magnetosphere is compensated by the same 

4 The black body radiation is assumed 



3. Discussion 

The basic conclusion of the investigation presented in this 
paper is that the observational appearance of a neutron 
star in the state of subsonic propeller strongly depends 
on the mode by which the envelope plasma enters the 
magnetosphere of the star. The canonical notion that the 
star in this state should be observed as a non-pulsating 
soft X-ray source of the luminosity i ssp is valid only in 
the particular case of the classical (Coulomb) diffusion. If, 
however, the penetration process is governed by the Bohm 
diffusion or/and the reconnections of the magnetic field 
lines, the mass accretion rate onto the star surface proves 
to be high enough for the accretion power to dominate 
the energy budget of the star. In this case the neutron 
star is expected to be observed as a steady low luminous 
(L x — L a ,diri + £ a ,rec) accretion-powered pulsar in which 
most of the X-ray emission comes from the hot (T cp ) spots 
of the radius ~ R cp situated in the magnetic pole regions 
on the surface of the neutron star (see Sect. 2.2). 

The entry of the accreting plasma into the interchange 
stable magnetosphere of a neutron star has been investi- 
gated by Eisner & Lamb ( |l984| ). They have shown that 
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there are no grounds for assuming that the Bohm diffu- 
sion or the field line reconnections modes of the plasma 
penetrating the magnetosphere are suppressed. They eval- 
uated the efficiency of the penetration process due to these 
modes to be of about 1% of the efficiency of that due to 
the interchange instabilities. This is too small to explain 
the observed luminosities of the brightest X-ray pulsars 
(~ 10 37 ~ 10 38 ergs _1 ) but is quite sufficient for the in- 
terpretation of the relatively low luminous X-ray sources 
(L K < 10 36 ergs^ 1 ). The efficiency of the plasma pene- 
trating the magnetosphere of a neutron star evaluated by 
Eisner & Lamb is comparable to that obtain ed in labora- 
tory experiments (e.g. Hamasaki et al. 1974 ) and derived 
from the observations of the solar wind enteri ng the mag- 
netosphere of the Earth (e.g. Gosling et al. 1991 ). The 
same results were obtained in the simulations of the so- 
lar wind penetration of the Earth's magnetic field due to 
diffusion (Winske & Omidi 1995 ) and reconnection of the 
field lines (Usadi et al. 1993 ). Thus the values £ = 0.1 
and a = 0.1 used in this paper seem quite reasonable. 
The estimation of the parameter A is more complicated 
and requires detailed study of the convective process in 
the envelope. This investigation, however, is beyond the 
scope of the present paper. Here I would like only to note 
that the estimate of the magnetic field scale by the size of 
a convection element is frequently used in the simulations 
of the magnetic field generation in astrophysical objects 
(e.g. accretion disks, active regions of the Sun, etc.). 

The last point I would like to address here is the de- 
termination of the state of a low luminous neutron star 
by the ratio between its magnetospheric and corotation 
radii, T = R m /R CO r- According to Eqs. @ and (JsJ) the 
value of T depends on the mass, the spin period and the 
dipole magnetic moment of the neutron star, and on the 
strength of the stellar wind in which the neutron star is 
situated, M c . The parameters of the neutron star can be 
evaluated if it is a component of a binary system and is 
observed as a pulsar with the cyclotron line feature in the 
X-ray spectrum. In this case the value of T can be ob- 
tained combining Eqs. (0), (§) and (|h: 



T = 



3£ 2/ V /7 



Ll /7 R 2 n i 7 (GM ns )^P? /3 ' 



(34) 



Here, £ = M a / M c is the efficiency of the accreting plasma 
penetrating the magnetic field of the star at the magne- 
tospheric boundary. 

As a rule, the value of £ is assumed to be equal to unity. 
This assumption is rather well justified if the neutron star 
undergoes disk accretion or/and if its X-ray luminosity 
exceeds the critical value, 



If, however, the plasma penetration into the magneto- 
sphere of a spherically accreting neutron star is governed 
by the diffusion or/and the field lines reconnection pro- 
cesses, the assumption £ ~ 1 is not valid. In particular, 
in the case of subsonic propeller £ Ss 0.01 and hence the 
value of T proves to be at least a factor of four smaller 
than that evaluated using the traditional method: 



T(£) = 0.27 



0.01 



2/7 



T(£ = l) 



This indicates that the overestimation of the parameter 
£ may lead to the mistaken classification of an accretion- 
powered source as supersonic propeller. An example of 
such a situation in the analysis of the low luminous state of 
the neutron star in the Be/X-ray transient A0535+26 has 
been recently presented by the author (Ikhsanov 20011] ). 
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3 10 36 a4q 4 «<' 



V 2 R-^eigs- 1 , 



at which the cooling of the plasma over the magneto- 
spheric boundary is effective and the boun dary is inter- 
change unstable (see Ghosh & Lamb 197S ; Arons & Lea 
1976a|). 



